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Abstract 

In this contribution we prove that the entropy of an N-body isolated system can not decrease and 
the entropy production should be non-negative provided the system possesses an equilibrium state. 
We define the entropy as a functional of the set of n-particle reduced density operators (n < N) 
generalizing the von Neumann fine-grained entropy formula. Additionally, as a consequence of our 
analysis we find the expression of the equilibrium n-particle reduced density operators which enter 
the definition of the entropy as well as the dissipated energy in an irreversible process. 
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Introduction.- It is a widely recognized fact that a general mathematical theoretical proof 
of the second law is still lacking. As stated in Ref. [l| and quoted here just for illustration 
sake "To the best of our knowledge no theoretical mathematical derivation of the second 
law has been given up until now; instead it has been based on Kelvin's or Clausius's prin- 
ciples of the impossibility of perpetual motion of the second kind P], which are based on 
experiment Isll " . This lack of definitive theoretical proof has lead to reports on the violation 

n n n 

of the second law |^ or tests over its validity in some particular casesp, p. 

The first significant contribution to the interpretation of the second law of Thermody- 
namics and the explanation of irreversibility goes back to Boltzmann. Nevertheless, it is 
known that Boltzmann's contribution was criticized by arguing that this contradicts the 
predictions based on the microscopic equations of motion. Later, Gibbs and P. Ehrenfest & 
T. Ehrenfest worked on this problem by introducing coarse-graining. However, those coarse- 
graining analyses require the introduction of a priori equal probability principles, which are 
hard to justify on physical grounds as was criticized by Einstein [3]. 

In this scenario, generalizing the Gibbs-von Neumann's statistics, it is our contention to 
provide a mathematical theoretical proof of the second law for an adiabatically isolated sys- 
tem based on first principles without coarse-graining. Our starting point is the description of 
the state of an isolated N-body system in terms of the set of n-particle reduced density opera- 
;ors in the framework of the BBGKY [Bogolyubov-Born-Green-Kirkwood-Yvon] description 
^. Notice that a global description of the system in terms of the full N-particle density op- 
erator is justified for equilibrium systems. Nonetheless, the nonuniformity of nonequilibrium 
systems leads to a random clusterization and the same global quantity may be obtained for 
an infinite number of realizations. This fact is taken into account in the BBGKY hierarchy 
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making this an appropriate framework for the description of nonequilibrium systems 
this context, since the collisions are made explicit through the collision term in the equations 
of motion, the n-particle reduced density operators are not constants of motion, therefore 
a way of defining the entropy to embody irreversibility might be in terms of this set of 
reduced density operators. This is what we do here: we propose a functional of the set of 
n-particle reduced density operators which generalizes the von Neumann relative entropy as 
the nonequilibrium entropy of the isolated N-body system. We will show that this entropy 
can not decrease and its rate of change or entropy production, should be non-negative. In 
addition, as a consequence of obtaining the entropy production we find the expression of 
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the equilibrium n-particle reduced density operators. In other words, we can know nothing 
about the equihbrium state without previously solving the dynamics of the system, an idea 
which coincides with Einstein's point of view on the subject 7|. 

Hamiltonian Dynamics. - Let's consider a dynamical system of N identical particles whose 
Hamiltonian is given by the N-particle Hamiltonian 

N ^ N 

j=i j^k=i 

where H(j) denotes the individual energy of the j-th particle and (j, k) the interaction 
energy between the j-th and /c-th particles. The state of the system is completely specified 
at a given time by the N-particle density operator D*^^)(l, N) which evolves according to 
the Liouville-von Neumann equation [q] 



^n^D%^ = L.,,,sD\T , (2) 



dt 

where -D^^'* = (a|D^^) are matrix elements and the Liouvillian operator L is defined 
through 

Lapas = [h^Is^^M - H^^^^6aS^ . (3) 

Here, repeated Greek indices mean summation. However, a strictly equivalent alternative 
description of the state of the system can be given in terms of the set of n-particle reduced 
density operators 3], 

I?={D(°),D«, ,D(^)} , (4) 

where T> constitutes a density superoperator whose components are the n-particle reduced 
density operators D*^")(l, ....,n) obtained by taking the trace over N — n particles, namely, 

and D^^p = 5af3 {i-e. D^^^ coincides with the unit matrix). These operators satisfy the 
normalization condition 

On the other hand, the dynamics of the D^"-' follows after performing the operation Tr 

{n+l,...,N) 

on both sides of Eq. Thus, we obtain the so-called quantum BBGKY hierarchy of 



equations 




(7) 



where H*^"^ is the n-particle Hamihonian 



(8) 



j=i j^k=i 



and 



(9) 



The second part on the left hand side of Eq. ((7j) constitutes the streaming term which gives 
the unitary evolution while the contribution on the right hand side is the collision term 
responsible for the non-unitary evolution of D^^^^ Streaming and collision terms balance out 
at equilibrium when D*^") = D^g''. 

Entropy and irreversibility.- In view of the fact that the statistical state of the system is 
given in terms of the set of reduced density operators, we define the nonequilibrium entropy 



a functional of the density superoperator V which constitutes a generalization of the relative 
or conditional von Neumann entropy. This entropy is the sum of all the n-particle entropies, 
a definition which is not redundant since the nonequilibrium system behaves as a random 
mixture of n-particle systems. Here Seq is the equilibrium entropy which corresponds to the 

equilibrium density superoperator V^q (=|d'^°), D^eq , , De^^j). This equilibrium density 

superoperator is obtined when the right hand side Eq. ((7j) is balanced by the second term on 
the left hand side of this equation, which occurs in a time scale related to the collisions 
Moreover, Seq must coincides with the thermodynamic entropy given by the Clausius's, 
Boltzmann's, Gibbs's etc. entropies. 

It should be pointed out that in defining the nonequilibrium entropy through Eq. , we 
assume that the BBGKY hierarchy describes a mixture of 'compressible' fluids, each consti- 
tuted by the set of n-particle clusters of the same size n. The origin of this 'compressibility' 
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keTT {V (InP- InVeq)} + S, 



eg 




(10) 



n=0 
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is the interaction between different ffuids which leads to the creation of n-particles clusters 
at the expense of the annihilation of p-particle clusters with n ^ p. Each fluid contributes 
with its entropy, the n-particle entropy, to the total nonequilibrium entropy of our system. 
Moreover, an additional consequence of the molecular collisions is the fact that each of the 
n-particle entropies can fluctuate and therefore also S. In the thermodynamic limit the 
fluctuations of S mentioned here should be gaussian, nevertheless this is not necessarily the 
case for finite N . 

More interestingly here the most important property of this entropy is its direction of 
change in a natural process. To elucidate this, we first note that since XlnX is a convex 



function and D*^") and Dig-* both are density operators 
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Tr {D(") (lnD(") -lnD(^))} > (11) 

which given D*^") = ^^Dm^ |m)^"''"-"-* (m| and oig^ = X]r-^eg,r jr)^"''^"'' (r| where d'^ and 
D^oi^r are eigenvalues and \m)^'^^ and |r) eigenvectors, can be written 

{^-"Hr I m)Wf (inZ^W _ inZ^W )| > . (12) 

Additionally since Dm^ gives the probability of the state \m) at time t, it should be bounded 
so that 

< 1^"^^ I [Dt^ (In/^L"^ -InZ^Si)} < ^ (13) 

q,m 

for a given constant C, which implies that the entropy S possesses a lower bound 

TV N 
n=0 q,m n=0 

Thus, one concludes that the nonequilibrium entropy can not decrease. 

Entropy production. - Following the previous line of thought, by using the BBGKY hierar- 
chy Eq. (0) one can compute the rate of change of S Eq. (|TO|l . i.e. the entropy production, 
obtaining 

^^ = _,riA;B$^^Tr{[H(")+^("),D(")] (inD^:;) - InD^"))} = 

CJT' 71/. 

n=0 

- ih-'kB Y ^Tr { [(H(") + ^(")) , InDg)] D^")} > (15) 

n=0 ^' 
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where the cydic invariance property of the trace has been used and [..,••] is the commutator. 
Moreover, JF^") has been defined through D(")jF(")= Tr |_D("+i)F*^"+^\ ) does accounting 

(n+l) 

for the interaction among the fixed particles of the n — th cluster with the remaining N — n 
particles of the system. The entropy production (jl5|) is zero at equilibrium when D*^"^ = D^^^ 
whereas in any other case it is generally non-zero, so in view of our conclusion after Eq. (|14p 
the entropy production should be positive. Thus, it might be said that this constitutes a 
general microscopic proof of the second law. 

In addition, as a consequence of Eq. (fT3j). one obtains that 

[H(")+^W,lnD(^)] =0 (16) 

which suffices to satisfy the extremum condition S{dS/dt) /SY)^"^^ \eq= meaning that 

lnDW = ir(H(")+^W) (17) 

for a given proportionality constant K, which after substitution in Eq. (fT^ allows us to 
rewrite this equation 

TV 



|5 = -^r^fc^i^ 5^ iTr { dW} , (18) 

n=0 

a compact way of writing dS/dt. Our result given through Eqs. (fTH|) and (fT^ is exact, 
thus possessing the generality that is lacking in the Boltzmann's H-theorem restricted to 
the molecular chaos regime. 

In the case that the system is in contact with a heat bath one can say that T>eq{t) and 
Seq{t), thus we can write 

where a is the irreversible entropy production given by right hand side of Eq. (jlSp and 
Q/To = dSeq/dt, with Q being the the rate of heat interchange and To is the bath temper- 
ature. According to our previous discussion a > 0, therefore 

6S>^^ (20) 

which corresponds to the Clausius' form of the second law. Hence, 

ToScT = ToSS - SQ (21) 
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coincides with the energy dissipated in the irreversible process. On the other hand, for 
reversible processes T>{t) = 'Deq{t) at every moment, therefore a = and Eq. pUj) becomes 

SS=^^ (22) 

Conclusions - We find that the description of the N-body system in the framework of the 
BBGKY hierarchy enables us to prove that the nonequilibrium entropy increases in a natural 
process, which constitutes a general microscopic proof of the second law. We emphasize that 
the entropy should be defined as a functional of the density superoperator. Likewise, our 
theoretical analysis allows us to obtain the expression of the set of equilibrium n-particle 
density operators, i.e. the equilibrium density superoperator. This fact shows that the state 
of equilibrium is determained as a consequence of the dynamics and not given a priori. 

As to how irreversibility can be derived from the reversible equations of mechanics, al- 
though the BBGKY hierarchy (jZj) expressing the microscopic dynamics is invariant under 
the time reversal operation, thus representing a set of reversible equations, this fact is not 
contradicted by the law of increase of nonequilibrium entropy we have derived. Let me 
explain, according to Eq. ((Tj), a stationary state is reached when the right hand side of this 
equation is balanced by the second term on the left hand side, this being the equilibrium 
state since this is the only stationary state for an . This stationarity or 

equilibrium condition is the same both for positive and negative time without distinction 
since we could say that at equilibrium there is no distinction between positive and negative 
time because at equilibrium time does not exist. Thus, the evolution along the time line 
either in the direct or inverse sense always ends in one and the same equilibrium state. 
Additionally, we have shown in the discussion after Eq. (|TO|l that the process of relaxation 
towards equilibrium entails the production of entropy according to Eqs. and ()18j) . thus 
having the character of an irreversible process. In conclusion, it could be said that time in- 
version is a mathematical artifice since there is only one Universe whose state given through 
the density superoperator (distribution functions in the classical language) evolves in one di- 
rection, the direction along which entropy increases up to its equilibrium value S^eq, wether 
this direction be direct or inverse. According to the time reversal invariance this is what 
one would predict. In fact, let us suppose that we have two identical isolated systems A and 
B, which at t = have positive and negative velocities, respectively . In those conditions, 
the past of B is the future of A or conversely, the future of B is the past of A. This can 
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be reworded in the following terms, if at a given time let us say t — 0, an isolated system 
is out of the equihbrium state, then 5" is at a local minimum and thus both dS/dt > for 
t > and dS/dt < ior t < 0, hence it follows that dS/dt is not necessarily a continuous 
function of time. Therefore, in those terms there is no reversible paradox. 

It should be stressed that in the basis of our treatment there is the assumption that 
the state of the system at a given time is determined by an statistical operator which is 
a function of the wave function, and therefore the Poincare recurrence theorem does not 
apply. In the classical language we would say that we do not deal with trajectories in the 
phase space but rather with distribution functions which determine the state of the system 
at a given time. A more fundamental proof of the fact that our theory does not incur in the 
recurrence paradox comes from the behavior of the entropy production dS/ dt as a function 
of the time explained above. 

To conclude, we have shown here that the reversible microscopic dynamics induces macro- 
scopic irreversibility. 
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